A glance at the imaginary world of ultracold atoms 
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From ultracold atoms to quantum chromodynamics, reliable ab initio studies of strongly interact- 
ing fermions require numerical methods, typically in some form of quantum Monte Carlo. Unfortu- 
nately, (non-)relativistic systems at finite density (spin polarization) generally have a sign problem, 
such that those ab initio calculations are impractical. It is well known, however, that in the rela- 
tivistic case imaginary chemical potentials solve this problem, assuming the data can be analytically 
continued to the real axis. Is this feasible for non-relativistic systems? Are the interesting features of 
the phase diagram accessible in this manner? Introducing complex chemical potentials, for real total 
particle number and imaginary polarization, the sign problem is avoided in the non-relativistic case. 
To give a first answer to the above questions, we perform a mean-field study of the finite-temperature 
phase diagram of spin-1/2 fermions with imaginary polarization. 



Ultracold Fermi gases provide an accessible and clean 
environment to study quantum many-body phenom- 
ena [TJ [3J, ranging from Bose- Einstein condensation 
(BEC) to Bardccn-Cooper-Schrieffer (BCS) superfluid- 
ity. In the dilute limit, where the range of the interac- 
tion is smaller than any other scale, a single parameter 
(/c F a s ) _1 , where a s is the s-wave scattering length and 
k F the Fermi momentum, describes the microscopic in- 
teractions completely. These are tuned by an external 
magnetic field in the presence of a Feshbach resonance. 

At large a s (in the crossover between BEC and BCS), 
these systems display universal properties. Here, the 
scale for all physical observables is set by solely by A; F 
(or, equivalently, the density), which is the only scale 
left in the problem. Thus, no obvious small expansion 
parameter exists in this limit, which represents a major 
challenge for theoretical many-body approaches [2]. Re- 
cently, experiments in this so-called unitary regime have 
achieved high precision [4] , potentially allowing to bench- 
mark and improve the available theoretical methods. 

Despite tremendous experimental and theoretical ad- 
vances, our understanding of ultracold Fermi gases at 
unitarity remains incomplete, most notably for the case 
of spin-imbalanced systems. For large-enough imbalance, 
one expects systems to undergo a phase transition from 
a BCS-type superfluid to a polarized normal gas. Such 
a transition has been observed in experiments at MIT 
and Rice university [7J and is in accordance with various 
theoretical studies [5] |S] . 

Apart from ultracold g better understanding 

of spin-imbalanced systems is also of great importance 
for other research fields. For example, almost all stable 
nuclei have more neutrons than protons and therefore 
fall into this class of systems. As a consequence, lat- 
tice MC calculations of nuclei [TU] therefore suffer from 



similar problems as ab initio lattice calculations of spin- 
imbalanced Fermi gases. 

We shall focus on the a s — > oo limit for a spin- 
imbalanced two-component Fermi gas at zero and finite 
temperature. Unlike previous studies (5] |5], however, 
we consider a complex- valued chemical potential (i. In 
ab initio MC calculations, this allows to circumvent the 
sign problem which spoils studies with spin imbalance. 
This approach parallels that of purely imaginary \i in 
relativistic quantum field theories, where it enables the 
analysis of the phase structure of QCD [TT] at finite 
density on the lattice. In that case, the existence of a 
critical end-point of the line of chiral transitions in the 
T-fi plane is still an open question, and the answer ap- 
pears to be out of reach for imaginary-/! approaches |12) . 
For non-relativistic fermions in the BEC-BCS crossover, 
on the other hand, the tri-critical end-point of the line 
of (second-order) superfluid transitions is known to ex- 
ist. Moreover, our present study suggests that this point 
might be accessible to lattice calculations by the use of 
complex- valued chemical potentials. 

In this first analysis we employ a mean-field approach, 
as discussed elsewhere for the case of real-valued \i (see 
e -g- 0), to study the phase diagram in the complex- 
valued case. Although this can only be viewed as a 
lowest-order approximation, it relies only on a single in- 
put parameter (e.g. kp) as is the case for the full eval- 
uation of the associated path-integral using, e.g., MC 
calculations. Thus, our statements do not suffer from a 
parameter ambiguity but only from an uncertainty asso- 
ciated with the underlying approximation. This can be 
understood on very general grounds from an analysis of 
the fixed-point structure of fermionic theories |13j . 

We begin by discussing a few general aspects of non- 
relativistic theories with complex-valued chemical poten- 
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tials. In general, the grand canonical partition function Z 
of non-relativistic fermions reads 



Z(T,fi,h) =Tr 



-,9(iTH3(iV t +^)-/t(# t -JV 4 )) 
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where T is the temperature and j3 = 1/T. We shall as- 
sume that the Hamiltonian H describes the dynamics of a 
theory with two fermion species, denoted by f and J,, in- 
teracting via a two-body interaction. The operators iV*. i 
denote the particle number associated with each species, 
and /j^^ are the corresponding chemical potentials. For 
convenience we have immediately introduced the average 
chemical potential /2 = (p^ + /Xj_)/2 and the asymmetry 
parameter h — (p^ — /i^)/2. 

As is well known, MC calculations for unitary fermions 
can be performed without a sign problem for h = (see 
e.g. Refs. 5, 6 ). This is not true in general, however, as 
polarization yields a sign problem regardless of the form 
of the interaction. To proceed, we consider an imaginary- 
valued asymmetry parameter h, which corresponds to 
studying a theory with complex- valued p^ ^ and there- 
fore define h = where /i : is a real quantity. It is easy 
to verify that MC calculations with imaginary-valued 
asymmetry can be studied with standard methods with- 
out a sign problem: the fermion determinants appearing 
in the probability measure are complex conjugates of one 
another. By analytically continuing Z(T, p, one ob- 
tains Z(T, p,h), which is the central quantity in studies 
of imbalanced Fermi gases. 

To understand whether the tri-critical end-point is ac- 
cessible with such an approach, we study the mean-field 
phase diagram with complex-valued chemical potentials. 
We compute the mean-field potential for the U(l) order- 
parameter, from the path-integral representation of Z: 
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where 
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and tp T — (ip^jip^) and g denotes the bare four- fermion 
coupling. The dimensionless renormalized four-fermion 
coupling g^gA is related to the scattering length a s by 



AnAg 1 = 



-rcg 
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Here A denotes the ultraviolet (UV) cutoff and the con- 
stant c rog . > depends on the regularization scheme. 
We use units such that 2m = 1, where m is the fermion 
mass. We represent the interaction via an auxiliary 
scalar field ip ~ <j> ^f^j.) where the parameter g v is cho- 
sen to reproduce the four-fermion term in the action. 
Fermions can then be integrated out straightforwardly, 
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Figure 1. (color online) Phase diagram of an ultracold Fermi 
gas at unitarity in the (T,h) plane. The solid (black) curve 
is a line of second-order phase transitions, which ends at a 
tri-critical point (h cp / p,T cp / p) and is followed by a line of 
first-order transitions (see e.g. Refs. [5]). The (red) dashed 
line is T/p — (T cp /h cp )h/ p and the (green) dashed-dotted 
line is TrT/p — h/p. The (light-blue) thin curves are possible 
analytic continuations obtained from Pade approximants of 
order -/V max = 1, 2, . . . , 5 (indicated by the numbers), see text. 
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Figure 2. (color online) Phase diagram in the {T,h{) plane. 
The solid line is a line of second-order phase transitions be- 
low which the fermion gap is finite. The (red) dashed and 
(green) dash-dotted lines are as in Fig. [I] with h T replacing h. 
Between the (red) dashed lines and the (green) dash-dotted 
lines there is a region with T c — > oo for p > 0. For p < 0, 
however, T c remains finite. 



as the resulting action is quadratic in those fields, which 
yields the order-parameter potential: 



pU(ip) = -2/3/2M : 



d 3 q 
(2^)3 
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This potential is directly related to the grand canonical 
potential: = VU(p ), where V is the volume of the 
system and <p denotes the mean-field value of ip that 
minimizes the potential. In that state, 5 2 |v? | 2 can be 
identified with the fermion gap A. The latter is the or- 
der parameter for spontaneous U(l) breaking associated 
with a supcrfluid ground state. From the (regularized) 
grand canonical potential we can derive all thermody- 
namic observables. Our results for dimensionless (uni- 
versal) quantities in the unitary limit, such as the crit- 
ical temperature T c /p for the superfluid transition, the 
gap A//1 2 or the ground-state energy E/p,, are indepen- 
dent of p and g v , as expected. To compute the critical 
temperature T c , it is convenient to employ the gap equa- 
tion (dU /dip)\ v = and exploit the fact that the fermion 
gap A ~ (pQ vanishes identically at T = T c . 

From Eq. Q, it is apparent that U is 27r-periodic in 
/3/ij when we replace h with i/q. This is not only true 
for the mean-field approximation but also for the full 
theory, as can be seen from Eq. ^ by taking into ac- 
count that h effectively shifts the Matsubara modes of 
the fermions v n = (2n + 1)ttT. [Note that, loosely speak- 
ing, d T is replaced by \v n when we formulate the action S 
in momentum space.] It then follows that it is not possi- 
ble to study arbitrary asymmetries with this technique: 
/q is bound to values /3/q < ir. However, the proposed 
approach still allows to study a large part of the phase 
diagram in the physical (T, ft,)-plane. We show next that 
our mean-field results suggest that a (tri)critical point in 
this phase diagram might indeed be accessible in lattice 
MC calculations with imaginary h. 

In Fig. [I] we show the well-known mean-field phase di- 
agram in the (T,h) plane (see e.g. Refs. [9 ). Rather 
than discussing the appearance of inhomogeneous phases 
(Sarma and/or FFLO), we only discuss the phase tran- 
sition lines for homogeneous phases. In Fig. [2] we show 
the phase diagram in the plane spanned by the tempera- 
ture and the imaginary-valued polarization. As discussed 
above, this phase diagram is 27r-periodic in /3/q. We 
therefore only show the domain /3/q G [— tt, it], which 
is indicated by the blue dashed lines both figures. As- 
suming one first computes the phase diagram in Fig. [2] 
e.g. via lattice MC, we conclude that we would then only 
have access to the temperature regime T > h/it in the 
physical (T, h) plane of Fig. [TJ Nevertheless, this repre- 
sents a fairly large part of the phase diagram, which is at 
the heart of theoretical and experimental studies. Most 
remarkably, our analysis suggests that the (tri) critical 
point lies within this regime of the phase diagram, i.e. 
that this point is within the reach of lattice MC calcu- 
lations with an imaginary asymmetry parameter. The 
phase transition line can then be obtained from an ana- 
lytic continuation of the results for T c (/q), as we mention 
below. 

The phase structure of the theory in the (T, hj) plane 
(Fig. [2]) is intriguing. It can be shown analytically 



that T c -!-oo for = (27V + l)?r with N eZ and p > 0. 
This is similar to relativistic fermion models, such as the 
Gross-Neveu model in (1 + l)d 14 , and again, we ex- 
pect this to hold beyond mean field. In fact, instead 
of the original Matsubara modes v n = (2n + 1)tvT, we 
see from the action |2| that the Matsubara modes of 
the fermions effectively assume the form v n = 2mrT 
for /3/q = (2N + 1)tt usually associated with bosonic de- 
grees of freedom. Thus, the fermions acquire a (thermal) 
zero mode in this case which tends to condense, inde- 
pendently of the actual value of the temperature. Con- 
trary to relativistic fermion models, however, we find nu- 
merically that already T c — > 00 for \(3h l \ > |(/37q)oo| = 
2.397 .... In other words, there is always a fermion con- 
densate for |(/3/ij)oo| < |/3 /iil < t and the U(l) sym- 
metry is not restored at high temperatures in this do- 
main. For \f3hj\ < |(/3/ii)oo|, the phase transition is of 
second order. The upper bound |(/3/q)oo| will take a dif- 
ferent value beyond the mean-field approximation. In 
our numerical studies we find that the value |/? C p^cp| 
associated with the (tri) critical point is slightly lower 
than \(/3hi)oa\j m what appears to be merely a coinci- 
dence. By studying the weak-coupling regime, one even 
finds that the difference between \ft C ph C p\ and |(/3/q)oo| 
is larger in that regime, at least in mean field [15]. Note 
that the absence of a (tri) critical point in the (T, hi) 
plane as well as the absence of U(l) restoration in the 
domain | (Ph-^oo | < |/3/q| < tt does not imply their ab- 
sence for the corresponding real-valued asymmetries. In 
analogy to relativistic fermion models [14], the analytic 
continuation of the phase boundary gives only the cor- 
rect behavior up to the (tri) critical point. From the ana- 
lytic continuation of the (full) order-parameter potential, 
however, the phase diagram can be mapped out in the 
physical (T, /i)-plane in the whole region where j3h < n, 
including the line of first-order transitions. 

As the grand canonical partition function Z, the order- 
parameter potential U is invariant under h — > —h. This 
allows us to expand Z (and other physical quantities) 
in powers of (/3h) 2 . In mean-field, we find that the ra- 
dius of convergence for the grand canonical potential 
is r = |/3/i|max = 7r for A = and /2 > 0, but r > tt 
in the case of a finite gap A. These observations facili- 
tate the analytic continuation from imaginary- to real- 
valued asymmetry parameter. When performing MC 
calculations of ultracold Fermi gases, one has now sev- 
eral options to do the analytic continuation. For ex- 
ample, one may fit the data of a given observable O 
for a given fixed temperature T — l//3 to the ansatz 

O = J2n=Q X ^o \^o^i) 2n 'i where are constants de- 
termined by the fit to the data. ./V max represents the 
truncation order (whose value depends on the amount 
of data available). Moreover, we have assumed that O 
has been made dimensionless with, e.g., a suitably chosen 
power of p,. From a simple analytic continuation of this 
polynomial, one then obtains the dependence of O on h. 
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Within our analytic study, we can check the feasibility of 
this procedure. For instance, we find that the pressure 
for < /3 h<l (with T = l//3 « p,/2) can be recov- 
ered from a fit to the imaginary- h data with iV max = 2. 
From the pressure one can then in principle compute the 
ground-state energy as function of h (or, equivalently, 
the so-called Bertsch parameter). However, the zero- 
temperature values of physical observables for finite h/jl 
are obviously not directly accessible within such an ap- 
proach. The Bertsch parameter at h = 0, on the other 
hand, is known from lattice MC calculations to approach 
rapidly its zero-temperature value below the superfluid 
transition [6]. We see a similar behavior in mean field, 
even at finite h/p,. Therefore it is conceivable that a re- 
liable estimate for the Bertsch parameter at T = and 
finite polarization can be extracted from lattice calcu- 
lations with various different temperatures and a subse- 
quent analytic continuation of the data. 

Instead of a simple Taylor expansion, one may employ 
more elaborate fit functions such as Pade approximants, 
also used in lattice QCD studies [TT] . In Fig. [I] for exam- 
ple, we have reconstructed the phase boundary at real- 
valued asymmetry by fitting the phase transition line in 
the {Phi, /3/i)-plane with the function 

c i + ££ m r «4i-cos(/^i)r 

1 + E7=i x ^[i-cos( W ' 

where -/V max again defines the truncation order. The a^'s, 
fej's and the constant C are determined by the fit. This 
ansatz already respects the 2-7r-periodicity of the theory in 
phj and can be generalized to observables other than T c . 
In Fig. [I] we show our results for the critical tempera- 
ture T c from such a fit for N max = 1, 2, . . . , 5, see Ref. [15] 
for details. Finally, we note that the fits may be even fur- 
ther optimized by choosing even more elaborate sets of 
basis functions [IB] . 

We have completely left out the discussion of inho- 
mogeneous phases and whether these phases also appear 
in the (T, hj) plane. While such a discussion is left to 
future work, we do not expect an inhomogeneous con- 
densate <p (x) ~ e lq °' x to show up in its well-known 
form, where the (center-of-mass) momentum <fo is de- 
termined by the difference in the chemical potentials of 
the spin-up and spin-down fermions. From a very sim- 
ple point of view, one may naively expect that the so- 
lutions ip (x) of the quantum equation of motion turns 
into f (x) ~ e~ qo ' x for complex- valued chemical poten- 
tials and might then no longer define the ground state. 

We have discussed the possibility of studying polarized 
Fermi gases with the aid of complex- valued chemical po- 
tentials. While the latter are not required in analytic 
studies, they are in MC calculations which would other- 
wise be spoiled by the sign problem. We have argued that 
the (tri)critical point is in principle within reach in this 
framework and the zero-temperature limit of observables 



might be indirectly accessible as well. This work there- 
fore suggests that, together with the experimental data 
at hand, future ab initio MC calculations with complex- 
valued chemical potentials have the capacity to push our 
understanding of collective many-body phenomena to a 
new level. Our present study marks the starting point 
and can already be used to guide these calculations. 
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